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Abstract. Given a collection of A'' solutions of the (3 + 1) vacuum Einstein 
constraint equations which are asymptotically Euclidean, we show how to con- 
struct a new solution of the constraints which is itself asymptotically Eu- 
clidean, and which contains specified sub-regions of each of the N given solu- 
tions. This generalizes earlier work which handled the time-symmetric case, 
thus providing a construction of large classes of initial data for the many body 
problem in general relativity. 



1. Introduction. 

An important problem in any theory of gravitation is the description of the 
motion of many-body systems. Various approximation schemes have been proposed 
to analyze this in general relativity, but no rigorous treatment has been provided 
thus far. A first step towards a solution of this question is to provide wide classes 
of initial data which solve the general relativistic constraint equations, and which 
are relevant to the problem at hand. There exists a rich family of initial data sets 
modeling isolated gravitational systems, but the non-linear nature of the constraints 
must be addressed when attempting to incorporate several such systems into a 
single one. In recent work [4j we have shown how this can be done by a gluing 
construction, under the restrictive hypothesis of time-symmetry. The aim of this 
work is to remove this restriction. 

Given a Riemannian metric g and a symmetric (0, 2)-tensor K on an oriented 
three-manifold M, the Einstein constraints map can be written in the for nE 

R{g) - K,,K'^ + {K\f 

The condition then for {g, K) to be the first and second fundamental forms of 
M embedded in a Ricci-flat space-time is that the vacuum constraint equations 
^{g,K) = (0,0) be satisfied. 

We recall that (g, K) on the exterior £^ of a ball in constitutes an asymp- 
totically Euclidean end (to order V) provided there are coordinates in which, for 
multi-indices |a| < £ + 1, |/3| < £, 

(1) \^"{9^, - '5.,)(x)| = 0(|xrl"l-i), \d^K,,{^)\ - 0(|x|-l^'-2), 

where d denotes the partial derivative operator. Note that throughout the rest 
of this work, we require that I > 2. We say that {M,g,K) is asymptotically Eu- 
clidean if M is the union of a compact set and a finite number of ends, all of which 
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are asymptotically Euclidean for {g, K) in the sense defined above. One read- 
ily verifies that every asymptotically Euclidean (AE) end possesses a well-defined 
energy- momentum vector (m,p). 

The starting point for constructing initial data for an N body system is the 
choice of the bodies. Each body is separately designated by the choice of an AE 
initial data set, and by the specification of a fixed interior region in each such AE 
solution. To perform this construction, we also specify a vacuum AE end disjoint 
from the interior region in each of the N data sets. We choose a collection of points 
which roughly locate the N bodies on a fiducial flat background. We then construct 
a new initial data set on a manifold obtained by excising a neighborhood of infinity 
in each chosen end, and then gluing these into a manifold Afcxt which is K."^ with N 
disjoint balls, centered around the chosen points, removed. In this way, the bodies 
are made to interact, by gluing the various ends into a fixed end. We now construct 
initial data on the resulting manifold which: i) is identical to the initial data from 
the N bodies away from Moxt, hence containing N chosen regions isometric to the 
specified interior regions of the bodies; ii) solves the vacuum constraints on Afcxt; 
iii) is identical to a space-like slice of a Kerr space-time sufficiently far from the 
bodies; and iv) has the centers of the bodies in a configuration which is a scaled 
version of the chosen configuration, where the scale factor can be chosen arbitrarily 
above a certain threshold. We emphasize that we preserve the original solutions 
away from a neighborhood of the gluing region. If the bodies are vacuum to begin 
with, we produce a solution to the vacuum constraints everywhere, but we can 
allow the initial data sets to be non-vacuum away from the chosen AE end. 

The construction actually produces a family of solutions depending on a param- 
eter e; roughly speaking, the bodies cannot be arbitrarily close together, but must 
be separated by a distance above a certain threshold. For each e, the distance 
between the bodies is on the order 0(e~^); the smaller e is taken, the further apart 
the bodies will be, and thus the weaker will be their initial interaction. We note 
that the energy-momentum four-vector of il/oxt in the resulting solution tends to 
the sum of the four- vectors of the bodies as e — > 0+ . 

We now state our main theorem. For clarity of exposition, we assume that the 
initial data sets are C°°-smooth, and refer the interested reader to [BUI] to formulate 
the case of a finite degree of regularity. 



Theorem 1. For each k — 1,...,N, let {Ek, g'^ , K'^) be a three-dimensional AE 
end which solves the vacuum constraints ^{g^,K*') = 0, with time-like energy- 
momentum four-vector (rnk,Pk)- LetUk C Ek be a pre- compact neighborhood of the 

N 

boundary dEk, and let M^xt = \ U ^k, where Bi, . . . , Bj^ are pairwise disjoint 

k=l 

closed balls in M.^ . There is an eo > so that for < e < eo, there is a solution 
{Mcxt, ge, Kt) of the vacuum constraint equations, with one AE end, containing 

N 

the disjoint union [j (J/^, 5*^, if '^), so that the distances between distinct Uk are 

k=l 

0(e^^). Near infinity {Mcxt,9t,K^) is isometric to a space-like slice of a Kerr 

N 

metric, with the ADM energy-momentum {m{g^),p^) satisfying rn{gg)—Y^ rrik < e 

k=l 

N 



ana 



E Pk 

fc=l 



< e. 
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The construction described here ahows us to glue together any finite number of 
asymptotically Euclidean ends which solve the vacuum constraint equations, and 
the construction is local near infinity in each end; i.e., any given compact subset of 
the end can be realized isometrically in the final metric g^. The local nature of the 
construction implies that we can allow the N original solutions to have multiple 
ends, and we can also allow nonzero matter fields supported outside a neighborhood 
of infinity in the chosen ends, as indicated below. 

Corollary 1.1. Let {AIk,g'',K'^), k — 1,...,N, be three-dimensional initial data 
sets with vacuum AE ends Ek C Mk of respective time-like ADM energy-momentum 
(mfc,pfe). Let Uk D Mk\Ek be chosen subdomains with EkHUk precompact. There 
is an eQ > so that for < e < sq, there is an initial data set (M, g^, K^^) which 

N 

contains a region U isometric to [J {Uk,g'^), for which {M \ U, ge, Ke) has one AE 

k=l 

end, with the same properties as those of {Mcyit,ge,Kt) as in Theorem]^ 

2. Preliminaries 

2.1. Kerr-Schild Coordinates. In order to establish certain estimates needed for 
the gluing carried out below, it will be convenient to use the explicit Kerr-Schild 
form of the Kerr metric (cf., e.g., [2]): 

2mr^ ^( m\ 

where 77 is the Minkowski metric (we are using the signature (— , +, +, +)), with 

1 z 

(2) Oudx'^ — dx^ — -35 7 [f{xdx + ydy) + a{xdy — ydx)] — —dz , 

+ r 

and where f is defined implicitly as the solution of the equation 

~4 ~2l 1 , 2 , 2 2\ 22 n 

r — r [X + y + z — a ) — a z = . 

We will continue to use the term Kerr-Schild coordinates for a coordinate system 
which has been obtained by a Lorentz transformation from the above. 

2.2. Global Charges. In establishing the main results, we use the fact that the 
Kerr space-times yield a family of initial data sets which admit coordinates (such as 
Kerr-Schild) with sufficient approximate parity symmetry that allows the definition 
of angular momentum J and centre of mass c, in addition to the four-momentum 
(m, p). Indeed the Kerr initial data sets we use satisfy the Regge-Teitelboim asymp- 
totic conditions, which say that in suitable AE coordinates the following estimates 
also hold: 

(3) 

= 0(|x|-|'''-3). 

The space of initial data satisfying ([3]) is known to be dense in the space of vacuum 
AE data [5]. 

Using such a coordinate system, we can compute the energy and linear and 
angular momenta using flux integrals at infinity {dae is Euclidean surface measure, 



a"(g.,(x)-.9,,(-x)) =0(|x|-l«l-2), 5^(x,,(x)+if,,(-x)) 
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V is the Euclidean outward normal, and r = |x|): 



{r=R} 



giii" ) 



{r=R} ^ 

Ji = ^ lim / ^(A> -XV9jfe)i^''^''rf'Te 

OTT R^oo J — ' 
{r=R} J'*-' 

{r=R} 

Note that in the last term, we can replace g in the center integrand by {g — (7euci)- 
Taken together, these give a set of ten Poincare charges associated to the end. We 
emphasize that we do not impose condition ^ on the initial data for the bodies; 
rather, we show in Proposition 13.11 that we can modify the given vacuum end of 
each body to a vacuum end in Kerr, preserving the data set away from the end. 

We recall the relation between the charge integrals and the constraints. Indeed, 
these charges arise from integrating the constraints against elements of the cokernel 
of the linearized constraint operator. By linearizing at the Minkowski data, we have 
'^{gEnci + h,K) = D^{h,K) + Q{h,K), where Q{h,K) = h*d'^h + dh*dh + dh*K, 
where denotes some metric contraction of the tensor product. In Euclidean 
coordinates at the Minkowski data, D^{h, K) = I '^{—'^{Kij^j — Kjj^i)), '^{hij^ij — 

j ij 

Thus for any vector and scalar pair {Y,N) which satisfies D^*{Y,N) = 
(0, 0), we have as a consequence of integration by parts that 
(4) 

/ {Y,N)■^gEucl+h,K)d^ie = B{R)-BiRo)+ f {Y,N)-Q{hJ<)dii, 

J J{Ra<r<R} 
{Ro<r<R.} 

where dfie is Euclidean volume measure, and where 

BiR) = [ (i-2)Y\K,,-K'Aj)+N{Y^{h,j,,-hu,j)-Y,iN,^h,,~N,hu)i^' da^ 

By letting F be a Euclidean Killing vector field, or letting be a constant or 
a coordinate function , we can easily relate B{R) to one of the above surface 
integrals defining the ADM energy-momenta. 

2.3. Hamiltonian Formulation of the Poincare Charges. It will be conve- 
nient to use the Hamiltonian formulation of the Poincare charges, as we review 
now (see Appendix E of [5] and references therein) . Let ^ be a three-dimensional 
spacelike hypersurface in a four-dimensional Lorentzian space-time {^,g). Sup- 
pose that ^ contains an open set with a time coordinate t (with range not 
necessarily equal to R), as well as a "radial" coordinate r e [i?, oo), leading to 
local coordinate systems (t, r, i;"^), with (z)^) providing local coordinates on a two- 
dimensional sphere. We further require that fl ^ = {t = 0}. Assume that 
the metric g^i, approaches the Minkowski metric rj^u as r tends to infinity. Set 
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— dx° Adx^ A dx'^ A dx'^ , and dSap — gf^-, ■, ■)■ The Hamiltonian analysis 

of vacuum general relativity in [3] leads to the formula 



for the Hamiltonian H{S^, g, X) associated to the flow of a vector field X, assumed 
to be a Killing vector field for the Minkowski metric, where 

SttU"^^ = J_ g^-,V.(|detg|.9-^[V'")^^ + ^d^g"''^VX^'a, 
V|dctg| 

with V the Levi-Civita connection for the Minkowski metric, and detg = det(gpcr). 
The integration over doc^ is taken, as usual, as a limit of integrals over spheres 
tending to infinity. We let Greek indices run from to 3, with x^ = t, and yll^'^l = 
^{A^^" — A'''^). We note that under enough approximate parity (such as holds for 
the Kerr data we consider), the above integrals converge [6]. 

If .y, viewed as a hypersurface in (^, g), has first and second fundamental forms 
{g, K) satisfying ^ and ([3]), then for the appropriate choice of X, the Hamiltonian 
yields the energy and momenta defined earlier. Indeed, if we let X = we 
get the energy-momentum one-form: H{S^,g,-^^) = p^, where po = m. With 
X = a^^gfr + ^^'gfrr, we get H{y,g,X) — mc', and with X = djix^ we have 
H{,9',g,X) = J'- 

We note that these Killing vectors correspond to Killing Initial Data, or KIDs, 
for the Minkowski metric [I . In Gaussian coordinates about the future-pointing 
unit normal is and for Y tangent to S^, N-^ -t- y is Killing for the Minkowski 
metric if and only if {D^)*{Y,N) — (0,0), where is the linearization of the 
constraints operator at the Minkowski data. We let U"^ (Y, N) correspond to X — 

In fact the surface integrals for the Hamiltonian at finite radii can be related 
to the constraints operator ^{g,K) = {Ti,p) by the following identity, where X = 
N-^ + y^gfr corresponds to the KID {Y,N), and g is a quadratic form in (gij — 
dij, ^i-, Kij) with coefficients uniformly bounded in terms of bounds on gij and 

f \]''^{Y,N)dSo.p = f W{Y,N)dSc,^ 

J {x'>=0,r=R} J {x°=Q,r=Ro} 

1 

(5) 



{Y'n +Np + q) dpg 
{xO=0,Ro<r<R} 



3. Constructing A^-Body Initial Data by Gluing 

3.1. Many-Kerr Initial Data sets. In this section we show how to attach a 
collection of N Kerr ends into one AE end while still solving the vacuum constraints. 
A special case of Theorem [5] below, for very special configurations (e.g., identical 
bodies placed symmetrically about a center) has earlier been established in [6l 
Section 8.9]. 

For < a < 6, let r(y, a, &) = {x e R'^ : a < |x — y| < 6} be a coordinate annulus 
of inner radius a and outer radius b centred at y. 



6 PIOTR T. CHRUSCIEL, JUSTIN CORVINO, AND JAMES ISENBERG 

Theorem 2 (Many-Kerr initial data sets). For each k — 1, . . . , N , let {Ek,g^ , K'^) 
be a Kerr asymptotically Euclidean three dimensional end, with Poincare charges 

Qk ■= {'mk,Pk,mkCk,3k) ■ 

Here the first entry is the ADM mass, the second is the ADM momentum, is the 
centre of mass, and Jfe is the total angular momentum. Suppose that the closures 

N 

of the annuli T(ck, 1,3) are pairwisc disjoint, and that ^ {mk,Pk) is time-like. 

There exists Cq > such that for each < e < cq there is a vacuum asymptotically 
flat initial data set {E,g^,K^) containing isometrically 

N 

y (^(e-lcfc,e-^2e-l),5^if'=) . 
fc=l 

N 

Proof. Let mr = to/c- Since mr is non-zero, we can translate the origin of 

k=l 

N 

coordinates so that the center of mass ^ rrikCk vanishes. Let -B(ro) C be a 

k=l 

Euclidean ball of radius ro = 5 + max{|ci|, . . . , |cjv|}, centred at the origin. 

We now construct a family of data on Qq := B{ro) \ ( IJ r(cfe, 0,1)1. Let x 

be a smooth nondecreasing function so that x(i) =0 for f < 9/4 and x{t) = 1 for 
t > 11/4. Define {g^,Kf) on fig as follows: 

• On each r(cA,, 1, 2), we let {ge, K^) be equal to {g^^k, Ke,k), where {gc,k,K^^k) 
are initial data for a Kerr metric in Kerr-Schild coordinates with global 
charges 

Qe,k = {^rnk,epk,emkCk,e'^3k) ■ 

• On each r(cfe, 2, 3), 

(5e(x),.^e(x)) = (1 -x(|x-Cfc|))(ffe,/c,-?^€,fc) +X(|x-Cfe|)(c/Eucl,0) . 

• On r(0, ro — 1, ro), we let 

(ffe(x),^e(x)) = (l-x(|x|-ro + 3))(5Eucl,0)+x(|x|-ro + 3)(ffe,ext(x),Jre,ext(x)) . 

where (ffcjextCx), -K'e,ext(x)) is a Kerr metric with global charge 
<3e,ext = (e(^rnfe + (5m),e(^pfe-|-5p),me,ext'5c,e^^Jfe + e^J^ , 



where 

e := {6m, Sp, 5c, 53)gQc . 
6 is a compact convex set which will be specified below. 

• On r(0,0,ro-l)\ U r(cfc,l,3), let {g„K,) = (ffEuci,0). 

fe=i 

In what follows we will use the fact, which follows directly from the above def- 
inition and from the Kerr-Schild form of the Kerr metric, that for any A: € N, we 
have 

(6) \\{9e,Ke) - (fl'Eucl,0)||c''(no) ^ ' 
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for some constant C — C{k). We thus have $(g£, K^) — 0{e), and $(ge, Ke) = in 

N 

a neighborhood of dH., where = B{ro) \ [J T{ck, 0, 2). 

k=l 

The goal is to modify {ge,K^) by adding a smooth deformation {5g^,5K^) sup- 
ported in so that '^{ge + 5ge, Kc + 5Ke) = 0. The proof proceeds in two stages, 
as we now describe. Recah that the hnearized constraints operator D<^ at the 
Minkowski data has a ten-dimensional cokernel ICq = ker(£)$*) spanned by the 
KIDs mentioned earlier (Section 12. 3p : this cokernel is an obstruction to solving 
the full system ^{ge + Sge,Ke + SK^) — for Sg^ and SK^. However, {ge,Ke) is 
close to the Minkowski data, and the constraints ^{g^,K^) vanish on a neighbor- 
hood of 9r2, and so in particular ^{g^^K^) belongs to the appropriate weighted 
spaces used in [6j|8]. Therefore, by applying Theorem 2 of [8] or Corollary 5.11 
of [6], we can at least solve the equation up to cokernel: we can find, for each 
9 G Q a.s above, a smooth deformation {6gf,S^) which is supported in il, sat- 
isfies \\{Sg^,,SK^J\\ (73 (o) < Ce (C is independent of 9 and e), and also satisfies 

^{(je + Sg^,Ke + 5K^) G C^o, where ^ is a smooth weight function that vanishes on 
9ri. For instance, in the notation of [6] , ^ = "(/i^ , for ^ a smooth function which near 
9r2 takes the form -0 = e^'*/'*, where d is a defining function for the boundary 917, 
and where s > 0. The collection of L^((i/Xe)-projections of $(gc -I- Sg^,K^ + SK^) 
onto a basis for /Cg defines a map 8 9 i— >■ M.^'^ which is continuous in 9. We show 
that this map vanishes for some 9, and that will complete the proof; a rescaling 
then yields Theorem [2l 

We now focus on computing these projections, which can be thought of as balance 
equations as we shall see. The primary tools for these computations are integration 
by parts together with the flux integrals ^ and Note that the difference 

\dfie — dfig^l is 0(e), and that the Taylor expansion around the Minkowski data 
(5Euc1,0) yields 

$a + 6gt, + 6K!) - + 6gt - 5euc1, Ke + SK^) + 0{^). 

As in Section 12.31 integrating this quantity against the basis of KIDs we get the 
boundary integrands of (|4]), integrated over dD,, plus O(e^) terms. Alternatively we 
can use an analogue of the integration-by-parts formula ([S]) to derive the balance 
equations in a form analogous to (8.7) of [B]. In either case, recall that on dVl 
our solution up to cokernel agrees with Kerr data. For such data, the differences 
between the limiting surface integrals (which give the charges) and surface integrals 
at finite radius are given by the divergence theorem as in ([5]). On the other hand, 
the Kerr data solves the vacuum constraints, so that the integrand over the annulus 
in ^ reduces to q. In the case of the data (5£,oxtj ^^^e.cxt) under consideration, this 
results in a O(e^) term; note that we use Kerr-Schild coordinates, in which the 
metric components take the form (77^^ + 0(TOe,cxt))- 

Keeping all this in mind, we can now compute the balance equations. Working 
first with 6(1) = (0, 1), the KID corresponding to the Minkowskian Killing vector 
^ (cf. Section [2^ . we obtain 



— (6(1), $(g, + 5gl + ))L.(n) 



an 



(7) 



e6m -f O(e^) . 



8 



PIOTR T. CHRUSCIEL, JUSTIN CORVINO, AND JAMES ISENBERG 



Next choosing e(i_|_j) = (i9i,0), for each i ~ 1,2,3, to be the KID corresponding to 
the Minkowskian KilUng vector ^ , we calculate the balance equation to be 

1 



T^{Hi+^)^'^i9e + 5glK,+5K'i))LHn) = e5f + 0{e'). 

IDTT 

Now let e(4_|_j), for each i = 1,2,3, be the KID corresponding to the Minkowskian 
Killing vector t^+a;'^;att = this corresponds to iY,N) = (0, a;'). The 
boundary integral around each is thus 



iaB(cfc,l) JdB{ck,l) 

[ U"'5(0, - cDdSc^p + [ W^{0, ci)dSa0 

JdB{ck,l) JdB{ckA) 



0+O(e2) 



4/ W^{Q,l)dS^0+O{e'') 
JaB{ck,i) 



£mfe+0(e2) 

= emfe4 + 0{e^) , 

where the first integral in the second line vanishes, up to quadratic terms in the 
metric, by definition of the centre of mass; we have also used linearity of the global 
charges with respect to the KIDs, and the fact that the ADM mass is the Poincare 
charge associated to the KID (0, 1). It follows that the associated balance equation 
reads 

(9) Y^(e(4+i), + Sg^,, K, + SK^,))L2(n) = m^.ext^c' - e ^ mfe4 +0(e^) • 



Finally, let di = and let 6(7+,), i = 1,2,3, be the KIDs corresponding to the 
Minkowskian Killing vector eijex^df, thus {Y,N) = {eijex^de,0). We then have 

/ U'"^{e^jex'de,0)dSc.p = [ U""^ {eije{x^ - 4 + 4)de,O)dSa0 

JaB{ck,l) JdB{ck,l) 



[ U"^{eije{x^ -cl)de,O)dSa0 

JdB{ck,l) 



dB(ck,l) 



2J«+0(c2) 



U"^{di,O)dSc,0+O{e'') 

dB(ck,l) 



epi+0(e^) 

We conclude that (recall that ^j^t angular momentum e'^^f.Jk + eS3) 
(10) -i-(e(7+i),$(5e + Sg^K, + dK'j)ma) = eST - e^(cfc x p^)^ + 0{e^). 
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We let 8 = (0, 0, 0, ^ Cfe x pfe) + Bq, where Bq is a closed ball around the origin 

fc 

chosen so that Qe.cxt has time-like four-momentum. We can invoke now the Brouwer 
fixed point theorem, in a way similar to the proof of Theorem 8.10 of [6 , to conclude 
that there exist eq small enough so that the right-hand-sides of the balance equations 
([7|l- (fT0l) can be all made to vanish for all < e < eo, which is the desired result. 

□ 

3.2. Reduction to Kerr Asymptotics. We recall the well-known result of [S1[S] 
which states that any AE vacuum end satisfying the Regge-Teitelboim condition ([3]) 
(with time-like ADM four-momentum) can be deformed, outside of a compact set, 
to a new vacuum initial data set such that the data agrees with that of a suitably 
chosen space-like slice of a Kerr space-time, outside of a compact set. We emphasize 
that this deformation can be performed to preserve any given pre-compact subset 
of the end; therefore we can apply this deformation to the data sets for each of our 
N bodies, preserving as large a compact subset of the original data as we like. We 
note for our main theorem that the four-momentum of the resulting Kerr can be 
made as close to the original four-momentum of the end as we like by performing 
the deformation at larger coordinate radii. 

Here we give a significant generalization of the above result, which removes the 
Regge-Teitelboim assumption ([3|) in the above gluing: 

Proposition 3.1. Let {g, K) be AE vacuum initial data on the exterior E of a ball 
in which has time-like ADM four-momentum (m, p). Let e > 0. For sufficiently 
large R, there is a vacuum initial data set {g,K) on E so that on Er\{\x.\ < R} we 
have {g,K) = (g,K), and so that on {|x| > 2R}, {g,K) is identical to data from 
a space-like slice of a suitably chosen Kerr space-time. If (m -|- (5to, p -|- ^p) is the 
four-momentum of{g,K), then \5m\ < e and \Sp\ < e. If, moreover, (0i holds, then 
also \6c\ < e and \63\ < e. 

Proof. The proof is a minor modification of that used in [6l[8] to prove the result 
for the case in which the condition ^ holds. In fact the primary modifications 
which are needed have been introduced in the proof of Theorem [21 

Given a vacuum AE end, we let : Ai ^ Afi he the scaling 0i?(x) = i?x, 
and let (5^, K^) = {R"^ (t)*i^g , R-^(t)*i^K). We note that since Q may not hold, the 
decay might not be good enough to make the centre of mass and angular momentum 
well-defined. However, we have 

^ / T.(^j'k-iKy,gf,)Yiiy'^dae = R-' J J^i^Jk ~ K'eg,k)Yi v^da^ 

{r=l} {r=R} 

(11, - o{'^y 

Similarly, 

(12) R J [^x^(g,^.,-5,,,^>J_^(.g,,^feV-5,y)]da^ 

{r=l} 

These two estimates follow from (|4]), together with the fact that {g,K) solve the 
vacuum constraints. Indeed the estimate Q{h,K) = 0{\x\^^) implies that 



^Compare [7j for smoothing arguments. 
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{Y,N) ■ Q{h,K) = 0(|xl-3) for {Y,N) = {Yi,x^) = 0(|x|). Now applying Q and 
the constraints, we obtain B{K) = B{Rq) + 0{\ogR). 

We let e = R~^, so that we can use some of the notation from the proof of 
Theorem[2] As above, we use a cutoff function to glue {g^,K^) to Kerr data with 

charges Qe.cxt — (^({m + 6m), e(p + (5p), e(m + S'm)Sc, eSJ^ on the annulus Ai, so 

that in a neighborhood of the inner boundary of Ai the data is identically {g^, K^), 
and in a neighborhood of the outer boundary the data is identical to the Kerr data. 
Again, we parametrize a family of such data with 9 :— (6m, 6p, 6c, 63) G = Bq, 
where Bq is a closed ball around the origin, chosen so that the four-momentum 
of Qe,cxt is time- like. Let {g^,K^) be the resulting glued data. We note that 

||.ge-ffEuci|lc«+i(Ai) + ll^e||c'^(Ai) = 0(e), that $(5,, i^e) = 0(e), and that 
vanishes in a neighborhood of dAi. As above, we solve the vacuum constraints up 
to cokernel, so that $(5^ -t- 6gl,k^ + 6K^) e C^o- 

We again analyze the projection of $(ge -I- 6g^,Ki + 6Kf) onto /Co, using the 
notation from the proof in the last section, along with ([TT]) and : 

^{e(,),^{g, + 6gl,k, + 6Kt))L^(A,) = t5m + 0{e') 

^{e^^^^^,^{g, + 6glk, + 6Kt))L-(A,) = e5f + 0{e^) 

■^{Hi+^)^^i9c + 6gl,k, + 6K'^,))L^(A,) = e[{m + 6m)6c' + o{l)] 

^{e(r+,),^{g, + 6gl,k, + 6Kl))L-(A,) = e(<5/' + o(l)). 

The proof now follows from the Brouwer fixed point theorem as before, followed by 
scaling back to the annulus A^i. □ 

3.3. Proof of Theorem [H 

Proof. We first apply Proposition 13.11 to deform the data on each end Ek to a 
new vacuum initial data set which agrees with the data on each Uk, and outside 
a compact set agrees with data from a suitably chosen space-like slice of a Kerr 
space-time. A rescaling of all the metrics then reduces the problem to one in which 
all of the initial data sets are Kerrian outside of a Kerr-Schild coordinate-ball of 
radius one. The result follows now by applying Theorem [2l □ 
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